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PROCEDE ET DISPOSITIF D'IMAGERIE RADIOGRAPHIQUE 
TRIDIMENSIONNELLE A FAIBLE DOSE D'IRRADIATION. 

(5j) Pour realiser une image en reconstruction tridimen- 
sionnelle a faibie dose d'irradiation du squelette d'un patient 
(P), on prend simuttanement deux images radiographiques 
perpendiculaires de cette partie du squelette, par balayage, 
on repere surchaque image des reperes de controle prede- 
termines de chaque vertebre, puis on reconstitue un modele 
numerique de la vertebre a partir des reperes de controle et 
d'un modele numerique generique predetermine de cette 
vertebre. 



POUR LA RECONSTITUTION 
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Procede et dispositif d'imagerie radiographique pour la 

reconstruction tridimensionnelle a faible dose 

d 'irradiation , 



5 La pr£sente invention est relative aux process et 

dispositif s d'imagerie radiographique pour la 

reconstruction tridimensionnelle a faible dose 
d T irradiation. 

Plus particulierement, 1' invention concerne un 

10 procede d'imagerie radiographique pour la reconstruction 
tridimensionnelle i faible dose d' irradiation, adapte pour 
calculer un modele a trois dimensions d'au moins un objet 
predetermine a observer dans un champ d f observation, ce 
procede comprenant les etapes suivantes : 

15 (a) prendre au moins deux images radiographiques a 

deux dimensions du champ d' observation, respectivement 
selon deux directions de prise de vue non paralldles, 

(b) reperer, sur chaque image radiographique, des 
reperes de controle appartenant audit objet a observer, 

20 (c) determiner une position geometrique de chaque 

repere de controle dans un referentiel ci trois dimensions, 

(d) calculer la forme a trois dimensions d'un 
modele reprfesentant ledit objet a partir d'un modele 
g£nerique predetermine correspondant audit objet, ce modele 

25 generique comportant des reperes qui correspondent aux 
rep§res de controle identifies sur les images 
radiographiques, le modele calcule etant obtenu par 
deformation du module g6nerique de fagon que ledit modele 
calcule suive une forme la plus proche possible d'une 

30 isometrie du modele g£n6rique tout en maintenant en 
coincidence les reperes du modele generique deforme avec 
les reperes de controle reconstruits a l'etape (c) . 

Des procedes de reconstruction en trois dimensions 
du type susmentionne ont ete divulgues notamment par Abdel- 

35 Aziz et al. ("Direct linear transformation from comparator 
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coordinates into object space coordinates in close range 
photogrammetry", Proc. ASP/UI Symp. Close Range 
Photogrammetry, Urbana, Illinois, USA, 1971) et Marzan 
("Rational design for close range photogrammetry" , PhD 
thesis, Department of Civil Engineering, University of 
Illinois, Urbana-Champaign, USA, 1976) . 

Dans ces procedes, tous les reperes de controle 
sont des points de controle stereo-correspondants et au 
cours de 1 ' etape (b) , on positionne ces reperes dans 
l'espace au moyen d'un algorithme dit "transformation 
lineaire directe" (DLT) , utilise notamment par Andre et al. 
("Optimized vertical stereo base radiographic setup for the 
clinical three-dimensional reconstruction of the human 
spine", J. Biomech., 27, pp 1023-1035, 1994), Par ailleurs, 
dans ces procedes connus, 1 r 6tape (d) est une etape de 
krigeage consistant en une interpolation/extrapolation du 
modele generique de l'objet £ observer, qui donne des 
positions estimees d'un grand nombre de reperes du modele a 
trois dimensions de l'objet a observer en fonction des 
coordonnees mesurees des reperes de controle stereo- 
correspondants et en fonction de la geometrie du modele 
generique* Cette etape de krigeage a ete decrite notamment 
par Trochu ("A contouring program based on dual kriging 
interpolation", Eng. Comput . 9, pp 160-177, 1993). 

Ces procedes connus presentent 1 1 avantage de 
permettre la realisation d'un modele tridimensionnel du ou 
des objets a observer, tout en permettant de reduire 
1' emission de rayonnements ionisants vers le champ 
d' observation par rapport a une information 
tridimensionnelle bas6e sur la reconstruction de coupes 
scanner telles que pratiquees dans les instruments actuels. 
Le modele tridimensionnel peut ensuite etre affiche sous 
differents angles de vue, par exemple sur un ecran 
d'ordinateur. 

Mais ces procedes souffrent d'un manque de 
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precision et sont mal adaptes pour examiner convenablement 
un champ d' observation 6tendu tel que par exemple 
l f ensemble de la colonne vertebrale d'un patient. 

Par ailleurs, MITTON et al. (Medical & Biological 
Engineering & Computing 2000, Vol. 38, p 133-139) ont 
decrit un procede de reconstruction tridimensionnelle qui 
utilise des reperes non stereo-correspondants en plus des 
reperes stereo-correspondants, pour reconstruire une image 
a trois dimensions d'une vertebre isolee, a partir de deux 
images radiographiques de cette vertebre prises 
successivement sous deux angles differents. 

Ce procede n'est toutefois pas adapts pour une 
utilisation medicale, ou il est necessaire d f avoir une plus 
grande facilite de mise en oeuvre, une meilleure precision 
notamment sur des champs d 1 observation etendus, et une plus 
grande rapidity de mise en oeuvre. 

La presente invention a notamment pour but de 
pallier ces inconvenients . 

A cet effet, selon 1' invention, un procede du genre 
en question est caracterise en ce qu'au cours de 1 1 etape 
(a) , les deux images radiographiques sont prises 
simultanement, par balayage, en deplagant en synchronisme, 
dans une meme direction de translation non parallele aux 
directions de prises de vues, au moins une source 
radioactive emettant deux faisceaux de rayons ionisants 
respectivement dans les deux directions de prise de vue. 

Le champ d 1 observation mentionne ci-dessus peut 
comprendre notamment le rachis, le bassin, ou encore le 
genou d'un patient, ou plus generalement etre constitue par 
tout ou partie du squelette du patient. Dans ces differents 
cas, les objets a observer peuvent etre constitues 
notamment par les os du patient compris dans le champ 
d'observation ainsi que la duree des operations de prise de 
vue. 

Grace aux dispositions susmentionnees, on obtient 
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une bonne precision de la reconstruction a trois 
dimensions, y compris pour des champs d' observation 
etendus, et ce en limitant la dose de rayonnements emise 
vers le champ d' observation . 

Cette precision est obtenue grace a la simultaneity 
des deux prises de vues, et grace a la prise de vue par 
balayage qui ameliore la precision dans la direction du 
balayage notamment pour les champs d' observation etendus. 

Dans des modes de realisation preferes du procede 
selon 1' invention, on peut eventuellement avoir recours en 
outre a l'une et/ou a 1' autre des dispositions suivantes : 

- au cours de l'etape (b) , certains des reperes de 
controle identifies sont des reperes de controle non 
stereo-correspondants visibles et identifies sur une seule 
image, et au cours de l'etape (c) , la position geometrique 
de chaque repere de controle non stereo-correspondant dans 
le referentiel a trois dimensions est estimee a partir du 
modele generique, en depla<?ant les reperes correspondant 
aux reperes de controle non stereo-correspondants du modele 
generique chacun sur une droite joignant : 

d'une part, la source radioactive a l'origine de 
l f image radiographique dans laquelle une projection de ce 
repere de controle non stereo-correspondant est visible et 
identifiable (la source radioactive etant positionnee sur 
1 1 emplacement de sa trajectoire qui correspond a la prise 
d* image ou est visible ledit repere de controle), 

et d' autre part, ladite projection de ce repere sur 
1' image radiographique, 

les reperes de controle non stereo-correspondants etant 
ainsi deplaces jusqu'a des positions respectives qui 
minimisent la deformation globale du modele generique de 
l'objet a observer / 

au cours de l'etape (c) , on minimise la valeur 
de la somme quadratique : 

1=1 
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ou X est un coefficient constant, m est un nombre entier 
de ressorts fictifs reliant chaque repere du modele 
generique A d'autres reperes de ce module, ki est une 
valeur de raideur predSterminee du ressort fictif d*indice 
if x io est la longueur du ressort fictif d'indice i dans le 
modele generique initial, et Xi et la longueur du ressort 
fictif d'indice i dans le modele generique deform^ 
(d'autres equations pour exprimer des deformations de 
lignes, surfaces ou volumes peuvent etre envisagees) ; 

au cours de l'etape (b) , certains des reperes de 
controle identifies sont des reperes de controle stereo- 
correspondants visibles et identifies sur les deux images, 
et au cours de 1'etape (c) , la position geometrique des 
reperes de controle stereo-correspondants est directement 
calculee a partir de mesures de position des projections de 
ces reperes sur les deux images ; 

- au cours de l'etape (b) , on repere sur chaque 
image radiographique des lignes de contour correspondant a 
des limites de l'objet observe et/ou a des lignes de plus 
grande densite optique a l'interieur desdites limites, ces 
lignes de contour comprenant des projections des reperes de 
controle sur les images radiographiques / 

- au cours de l T £tape (c), on determine des 
reperes du modele generique correspondant aux reperes de 
controle, lesdits reperes du modele generique comprenant 
des portions dudit modele generique qui se presentent 
tangentiellement par rapport aux rayons issus des sources 
radioactives et ayant genere les images radiographiques ; 

l f etape (c) comporte les sous-etapes suivantes : 
(cl) creer un modele generique recall en adaptant 
la taille du modele generique et la position de ce modele 
generique dans le referentiel, pour que les projections 
respectives du modele generique recale a partir des deux 
sources radioactives correspondent sensiblement aux deux 
images radiographiques, 
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(c2) selectionner des reperes du modele generique 
dont les projections sur au moins une des images 
radiographiques a partir de la source radioactive 
correspondante, sont les plus proches des lignes de contour 
reperees au cours de l'etape (b) , 

(c3) definir une surface enveloppe formee par des 
rayons issus de chaque source radioactive et ayant 
contribu£ a generer lesdites lignes de contour des images 
radiographiques, 

(c4) determiner certains reperes du modele 
generique recale correspondant a des surfaces dudit modele 
generique recale, qui sont tangentes auxdites surfaces 
enveloppes, les reperes du modele generique recale ainsi 
determines correspondant aux reperes de controle, 

(c5) et determiner la position geometrique de 
chaque repere de controle par projection du repere 
correspondant du modele generique recale sur la surface 
enveloppe correspondante ; 

les deux directions de prise de vue sont 
perpendiculaires l'une a I 1 autre ; 

- on fait emettre par chacune des sources 
radioactives un faisceau de rayonnements ionisants dans un 
plan perpendiculaire a la direction de translation ; 

les deux faisceaux de rayons ionisants sont emis 
respectivement par deux sources radioactives. 

Par ailleurs, 1' invention a egalement pour objet un 
dispositif d'imagerie radiographic pour la mise en oeuvre 
d'un precede tel que d6fini ci-dessus, ce dispositif 
comportant : 

des moyens d' emission de rayons ionisants 
comprenant au moins une source radioactive, ces moyens 
d ? emission etant adaptes pour emettre respectivement deux 
faisceaux de rayonnements ionisants vers un champ 
d' observation contenant au moins un objet a observer, dans 
deux directions de prise de vue non paralleles, lesdits 
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moyens d* emission etant deplagables simultanement selon une 
direction de translation non parallele aux directions de 
prise de vue, 

au moins deux dispositifs de detection disposes 
5 respectivement face aux deux faisceaux de rayons ionisants, 
au-dela du champ d' observation, pour mesurer les 
rayonnements ionisants ayant traverse ledit champ 
d' observation, ces deux dispositifs de detection etant 
deplagables en synchronisme avec les moyens d f emission dans 

10 ladite direction de translation, 

des moyens pour prendre simultanement au moins 
deux images radiographiques a deux dimensions du champ 
d 1 observation, par balayage simultan6 du champ 
d ? observation avec les sources radioactives et les 

15 detecteurs dans la direction de translation, 

des moyens d f identification pour identifier, sur 
chaque image radiographique, des reperes de contrdle 
predetermines appartenant audit objet a observer, 

des premiers moyens de reconstruction pour 

20 determiner une position geometrique de chaque rep£re de 
controle dans un r£ferentiel a trois dimensions, a partir 
d'un modele generique predetermine correspondant audit 
objet, ce modele generique comportant des reperes qui 
correspondent aux reperes de controle identifies sur les 

25 images radiographiques, 

et des seconds moyens de reconstruction pour 
calculer la forme a trois dimensions d f un modele 
repr^sentant ledit objet a partir du modele generique, 
lesdits seconds moyens de reconstruction 6tant adaptes pour 

30 determiner le modele calcule par deformation du modele 
generique de facon que ledit modele calcule suive une forme 
la plus proche possible d'une isometrie du modele generique 
tout en maintenant en coincidence les reperes du modele 
generique deforme avec les reperes de controle reconstruits 

35 par les premiers moyens de reconstruction. 
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Dans des modes de realisation preferes du 
dispositif selon 1 • invention, on peut eventuellement avoir 
recours en outre a l'une et/ou a 1 T autre des dispositions 
suivantes : 

les moyens d 1 identification sont adaptes pour 
identifier au moins certains reperes de controle visibles 
et identif iables sur une seule image radiographique dits 
repdres de controle non stereo-correspondants et les 
premiers moyens de reconstruction sont adaptes pour estimer 
la position g£ometrique des reperes de controle non stereo- 
correspondants en deplagant les reperes du modele generique 
correspondant aux reperes de controle non stereo- 
correspondants, chacun sur une droite joignant : 

d'une part, la source radioactive a l f origine de 
1' image radiographique ou une projection de ce repere de 
controle non stereo-correspondant est visible et 
identifiable, 

et d f autre part, ladite projection de ce repere sur 
l 1 image radiographique, 

les premiers moyens de reconstruction etant adaptes pour 
deplacer ainsi les reperes de controle non stereo- 
correspondant s jusqu'a des positions respectives qui 
minimisent la deformation globale du modele generique de 
I'objet a observer ; 

les moyens d 1 identif ication sont adaptes pour 
identifier sur les deux images des reperes de controle qui 
sont visibles et identif iables sur lesdites deux images, 
dits reperes de controle stereo-correspondants, et les 
premiers moyens de reconstruction sont adaptes pour 
determiner la position geometrique des reperes de controle 
stereo-correspondants par calcul a partir de mesures de 
position des projections de ces reperes sur les deux 
images ; 

chaque detecteur comprend une ligne de cellules 
de detection perpendiculaire a la direction de translation, 
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les faisceaux de rayonnements ionisants etant 
perpendiculaires a ladite la direction de translation ; 

les moyens d' emission et les detecteurs sont 
deplagables sur une distance de balayage d f au moins 70 cm ; 

les moyens d' emission comprennent deux sources 
radioactives a l'origine respectivement des deux faisceaux 
de rayons ionisants. 

D'autres caracteristiques et avantages de 
1" invention apparaitront au cours de la description 
suivante d'une de ses formes de realisation, donnee a titre 
d'exemple non limitatif, en regard des dessins joints. 
Sur les dessins : 

la figure 1 est une vue schematique d'un 
appareil de radiographic selon une forme de realisation de 
1' invention, permettant d'effectuer simultanement une prise 
de vue de face et une prise de vue de profil du patient, 

la figure 2 est une vue schematique en 
perspective d'une vertebre d'un patient examine au moyen de 
1' appareil de la figure 1, 

les figures 3 et 4 sont respectivement des vues 
de profil et de face de la vertebre de la figure 2, 
schematisant une partie des vues de profil et de face 
obtenues au moyen de 1' appareil de la figure 1, 

et la figure 5 est une vue en perspective 
representant un modele & trois dimensions de la colonne 
vertebrale et du bassin du patient examine au moyen de 
l f appareil de la figure 1, ce modele etant calcule a partir 
des vues de profil et de face prises au moyen de 1' appareil 
de la figure 1. 

Sur les differentes figures, les memes references 
designent des elements identiques ou similaires. 

La figure 1 represente un dispositif radiographique 
1 pour la reconstruction tridimensionnelle, comportant un 
bati mobile 2 deplagable verticalement de fagon motorisee 
sur des guides verticaux 3, dans une direction de 
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translation 3a. 

Ce bati entoure un champ d' observation 4 dans 
lequel peut prendre place un patient P debout. On peut 
ainsi observer la position des os du squelette de ce 
patient en station debout, ce qui est essentiel notamment 
pour les patients atteints de scoliose. 

Le bati mobile 2 porte une premiere source 
radioactive 5 et un premier detecteur 6 qui est dispose 
face a la source 5 au-dela du champ 4, et qui comporte une 
ligne horizontale 6a de cellules de detection. Le detecteur 
6 peut par exemple etre un detecteur gazeux sensible aux 
basses doses de rayonnements, par exemple tel que decrit 
dans le document FR-A-2 749 402 ou FR-A-2 754 068. Bien 
entendu, d'autres types de detecteurs, gazeux ou non, 
pourraient eventuellement etre utilises dans le cadre de la 
presente invention . 

La source radioactive 5 est adaptee pour emettre 
des rayons ionisants, notamment des rayons X, dans une 
direction de prise de vue 7 antero-posterieure par rapport 
au patient P, en traversant une fente horizontale 8 menagee 
dans un reticule 9 tel qu'une plaque metallique, pour 
g£nerer un faisceau horizontal 10 de rayonnements ionisants 
dans le champ d * observation 4. 

Par ailleurs, le bati mobile 2 porte egalement une 
deuxieme source radioactive 11 similaire a la source 5 et 
un deuxieme detecteur 12 similaire au detecteur 6, qui est 
dispose face a la source 11 au-delA du champ 4, et qui 
comporte une ligne horizontale 12a de cellules de 
detection. 

La source radioactive 11 est adaptee pour emettre 
des rayons ionisants, dans une direction de prise de vue 
laterale 13 par rapport au patient P, en traversant une 
fente horizontale 14 menagee dans un reticule 15 tel qu'une 
plaque metallique, pour gendrer un faisceau horizontal 16 
de rayonnements ionisants dans le champ d* observation 4. 
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On notera que les sources radioactives et les 
detecteurs pourraient le cas echeant etre en nombre 
superieur a 2, et que les directions de prises de vue de 
ces differentes sources radioactives pourraient le cas 
echeant ne pas etre perpendiculaires entre elles ni 
horizontales . 

Les deux detecteurs 6, 12 sont relies a un micro- 
ordinateur 17 ou autre systeme electronique de commande, 
equipe : 

d'une interface d' entree comprenant au raoins un 
clavier et generalement une souris (non representee), 

et d'une interface de sortie comprenant au moins 
un ecran 19 et generalement une imprimante (non 
representee) . 

Le micro-ordinateur 17 peut egalement etre relie 
aux moyens d' entrainement motorises (non representes) 
contenus dans les guides 3 et aux sources 5, 11, de fagon a 
commander le deplacement vertical du bati 2 et l f emission 
des rayonnements ionisants. 

Le dispositif qui vient d'etre decrit fonctionne 
comme suit. 

Au moyen du micro-ordinateur 17, on fait d' abord 
prendre deux images radiographiques du patient P, en 
faisant balayer le champ d'observation 4 par les faisceaux 
10, 16 de rayonnements ionisants sur la hauteur 
correspondant a la zone du patient a observer, par exemple 
le rachis et le bassin, voire 1' ensemble du squelette (a 
cet effet, le bati est de preference deplagable sur une 
hauteur d'au moins 70 cm, voire superieure aim). 

Au cours de ce mouvement, on enregistre dans la 
memoire du micro-ordinateur 17 deux images radiographiques 
numeriques, par exemple respectivement antero-posterieure 
et laterale de la partie examinee du patient, lesquelles 
images peuvent etre visualisees sur 1 1 ecran 19 du micro- 
ordinateur. 
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Chacune de ces images comprend gen6ralement 
plusieurs objets predetermines a examiner, par exemple des 
vertebres 20 telles que celle representee sch£matiquement 
sur la figure 2. 

5 Pour chacun de ces objets a examiner/ le micro- 

ordinateur 17 a en m6moire un modele generique a trois 
dimensions qui correspond a une forme moyenne de 1' objet en 
question, lequel modele generique est elabore a l'avance 
par des methodes statistiques en analysant un grand nombre 
10 d f objets similaires. 

Lors de l'affichage des images radiographiques sur 
1 ' ecran 19 du micro-ordinateur 17, le praticien peut par 
exemple indiquer au micro-ordinateur, notamment au moyen du 
clavier 18 ou de la souris, le type de chaque objet a 
15 examiner visible sur lesdites images, de fagon que le 
micro-ordinateur 17 determine le modele generique 
correspondant a cet objet. 

Par ailleurs, les modeles generiques utilises 
pourraient egalement etre constitues par des modeles 
20 prealablement realises par imagerie medicale sur le 
patient : dans ce cas, le procede selon 1' invention peut 
permettre par exemple de suivre Involution ulterieure du 
patient par des moyens plus simples, moins couteux et 
emettant moins de radiations que les moyens d' imagerie 
25 tridimensionnelle classiques. 

Le modele generique de chaque objet, par exemple de 
chaque vertebre 20 d'un squelette humain, comprend : 

les coordonnees d'une pluralite de reperes de 
controle, notamment des points correspondant a des reperes 
30 singuliers de cette vertebre, 

et les coordonnees d ? un grand nombre d'autres 
reperes de l 1 objet en question, par exemple au nombre 
d 1 environ 200 ou plus. 

Ces coordonnees peuvent etre exprimees par exemple 
35 dans un referentiel local X, Y, Z. Dans l 1 exemple 
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considere, l'axe Z correspond a la direction "axiale" de la 
colonne vertebrale, l'axe X est determine de fagon a 
definir avec l'axe Z le plan antero-posterieur de la 
vertebre 20, l'axe Y etant perpendiculaire aux axes X, Z 
susmentionnes. De plus, l'origine 0 du referentiel X, Y, Z 
est disposee au milieu des deux faces d' extremite axiales 
de la partie principale "tubulaire" de la vertebre, 
l f origine O etant par ailleurs positionnee pour que l'axe Z 
traverse la face axiale superieure de la partie principale 
de la vertebre en un repere CI tel que la distance de ce 
repere CI & 1' extremite avant C7 de ladite face axiale soit 
egale a environ 2/3 de la distance totale entre les 
extr<§mites avant C7 et arriere C8 de la section antero- 
posterieure de ladite face axiale superieure. 

Les differents reperes de controle C1-C25 
susmentionnes se repartissent en deux categories : 

- des reperes de controle "stereo-correspondants" 
C1-C6, visibles et identif iables a la fois sur 1' image 
radiographique laterale et sur 1' image antero-posterieure, 
ces reperes etant au nombre de 6 dans I'exemple considere 
(voir figures 3 et 4), 

et des reperes de controle "non stereo- 
correspondants" C7-C25, visibles et identif iables sur une 
seule image, ces reperes etant au nombre de 19 dans 
I'exemple considere, 

Le praticien identifie ces differents reperes de 
controle pour chaque objet a examiner (par exemple les 
vertebres et le bassin) sur chaque image radiographique, 
par exemple en "marquant" ces reperes a 1 1 ecran 19 par 
selection au moyen de la souris et/ou du clavier. De plus, 
les deux images sont calibrees, de fagon a pouvoir mesurer 
precisement la position de chaque repere de ces images dans 
un referentiel commun. 

Ensuite, on determine une position geometrique de 
chaque repere de controle de chaque objet, dans un 
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referentiel a trois dimensions, par exemple le referentiel 
X, Y, Z susmentionne ou un referentiel commun a 1' ensemble 
des objets a examiner. 

La position des reperes de controle stereo- 
5 correspondants C1-C6 est directement calculee a partir de 
la mesure de la position de ces points sur les deux images, 

De plus, la position geometrique de chaque repere 
de controle non stereo-correspondant C7-C25 dans le 
referentiel a trois dimensions est estimee a partir du 
10 modele generique, en deplagant chaque repere de controle 
st6reo-correspondant C1-C6 du modele generique jusqu'a sa 
position mesuree, et en d§plagant les reperes de controle 
non ster£o-correspondants C7-C25 du modele generique, 
chacun sur une droite joignant : 
15 . d'une part, la source radioactive 5, 6 a l'origine 

de 1' image radiographique ou une projection de ce repere de 
controle non stereo-correspondant est visible et 
identifiable, 

et d' autre part, ladite projection de ce repere sur 

20 1' image radiographique, 

les reperes de controle non stereo-correspondants etant 
ainsi deplaces jusqu'a des positions respectives qui 
minimisent la deformation globale du modele generique de 
l'objet a observer. 

25 En particulier, on peut minimiser ladite 

deformation en minimisant (par exemple au moyen d'une 
methode de gradient) la valeur de la sorame quadratique : 

m 

ou plus generalement S = X ^k r (x i - x i0 ) 2 f 

i=\ 

30 ou X est un coefficient constant predetermine, m est un 
nombre entier non nul representant un nombre de ressorts 
fictifs qui relient chaque repere de controle du module 
generique k d'autres reperes de controle, ki est un 
coefficient de raideur predetermine du ressort fictif 
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d'indice i, x io est la longueur du ressort fictif d'indice 
i dans le modele generique non d6forme, et Xi est la 
longueur du ressort fictif d'indice i dans le modele 
generique deforme* 

Enfin, on calcule la forme a trois dimensions d'un 
modele representant la vertebre 20 du patient, le modele 
calcule etant obtenu par deformation du module generique de 
fagon a maintenir la coincidence des points de controle du 
modele generique . deforme avec la position spatiale 
precedemment determine des points de controle et de fagon 
que ledit modele calcule suive une forme la plus proche 
possible d'une isometrie du modele generique, en 
travaillant cette fois sur 1' ensemble des points de modele 
generique . 

En particulier, l'obtention du modele a trois 
dimensions de chaque objet a examiner peut etre obtenu par 
le procede connu de krigeage ( 11 kriging") ♦ 

Apres le calcul du modele a trois dimensions des 
differents objets a examiner, le micro-ordinateur 17 peut 
assembler la totalite des modeles a trois dimensions des 
differents objets a examiner, en fonction de la position de 
ces differents modeles dans un referentiel absolu commun a 
tous ces objets, de fagon a obtenir un modele a trois 
dimensions comprenant par exemple 1 ? ensemble du rachis 21 
du patient et le bassin 22 de ce patient, comme repr£sente 
sur la figure 5. 

Une fois elabore, ce modele a trois dimensions peut 
etre presente sur l'ecran 19 du micro-ordinateur, ou 
imprime, sous 1' angle de vision voulu. Ce module d' ensemble 
peut egalement etre mis en mouvement a l'ecran en fonction 
des commandes du praticien. 

Le praticien dispose ainsi . d'un outil efficace 
d'examen pouvant servir a l f imagerie de toute partie 
notamment osseuse ou cartilagineuse du corps humain ou 
animal, et utile notamment pour le diagnostic des scolioses 
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ou pour le suivi pre ou post-operatoire lors 
d 1 interventions chirurgicales - 

Bien entendu, on peut 6galement calculer certains 
indices cliniques predetermines lies soit a la geometrie de 
5 l f ensemble examine, soit le cas echeant a la composition ou 
la densite des objets a examiner, estimees a partir des 
images radiographiques (cas de 1 1 osteoporose par exemple) . 

On notera que le dispositif radiographique 1 
pourrait le cas echeant §tre adapte pour 1 1 examen d T un 

10 patient couche, ce qui peut s'averer indispensable dans le 
domaine de la traumatologic . Dans ce cas, le patient P 
serait couche sur une table support, les faisceaux de 
rayons ionisants 10, 16 seraient chacun dans un plan 
vertical, et les sources 5, 11 se deplaceraient 

15 horizontalement avec les detecteurs 6, 12. 

Par ailleurs, il va de soi que dans tous les cas, 
le dispositif radiographique 1 peut etre utilise egalement 
en radiographic a deux dimensions, en plus de son 
utilisation en imagerie tridimensionnelle . 

20 On notera que le dispositif selon I 1 invention 

pourrait le cas echeant etre utilise dans des applications 
de radiologie non medicale. 

Par ailleurs, au lieu d'utiliser des reperes de 
controle C1-C25 definis a l'avance sur chaque modele 

25 generique, il serait possible de determiner et de 
positionner dans l'espace les reperes de controle a partir 
de lignes de contour de l'objet a observer visibles sur 
l'une ou 1' autre des deux images radiographiques. 

A cet effet, on pourrait en particulier proceder 

30 comme suit : 

on repere sur chaque image radiographique des 
lignes de contour correspondant a des limites de l'objet 
observe et/ou a des lignes de plus grande densite optique a 
l'interieur desdites limites, 

35 - on cree un modele generique recall en adaptant 
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la taille du modele generique et la position de ce modele 
generique dans le referentiel X, Y, Z pour que les 
projections respectives du module generique recale a partir 
des deux sources radioactives 5, 11 correspondent 
sensiblement aux deux images radiographiques, 

- on selectionne des reperes du modele generique 
dont les projections sur au moins une des images 
radiographiques a partir de la source radioactive 
correspondante, sont les plus proches des lignes de contour 
reperees au cours de l'etape (b) , 

- on definit une surface enveloppe formee par des 
rayons issus de chaque source radioactives et ayant 
contribue a generer lesdites lignes de contour des images 
radiographiques , 

on determine certains reperes du modele 
generique recale correspondant a des surfaces dudit modele 
generique recale, qui sont tangentes auxdites surfaces 
enveloppes, les reperes du modele generique recale ainsi 
determines correspondant aux reperes de controle, 

- on determine la position geometrique de chaque 
repere de contr61e par projection dudit repere de controle 
sur la surface enveloppe correspondante, 

puis on procede par exemple corame decrit 
precedemment pour reconstituer un modele a trois dimensions 
de l 1 ensemble de l'objet a observer, notamment par 
krigeage. 

On notera enfin qu'il serait possible de generer 
deux faisceaux ionisants non paralleles au moyen de deux 
reticules (par exemple deux fentes distinctes meraorisees 
dans une meme plaque metallique) a partir d'une source 
radioactive unique pour mettre en ceuvre 1' invention, en 
utilisant comme precedemment deux detecteurs disposes face 
aux deux faisceaux et deplagables en synchronisme avec la 
source et les reticules. 
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REVINDICATIONS 

1. Procede d f imager ie radiographique pour la 
reconstruction tridimensionnelle a faible dose 
d 1 irradiation, adapte pour calculer un modele a trois 
dimensions d ! au moins un objet predetermine (20) a observer 
dans un champ d' observation (4), ce procede comprenant les 
etapes suivantes : 

(a) prendre au moins deux images radiographiques a 
deux dimensions du champ d 1 observation, respectivement 
selon deux directions de prise de vue (7, 13) non 
paralleles, 

(b) reperer, sur chaque image radiographique, des 
reperes de controle (C1-C25) appartenant audit objet a 
observer, 

(c) determiner une position geometrique de chaque 
repere de controle, dans un referentiel a trois dimensions, 

(d) calculer la forme a trois dimensions d'un 
modele representant ledit objet a partir d f un modele 
generique predetermine correspondant audit objet, ce modele 
generique comportant des reperes qui correspondent aux 
reperes de controle identifies sur les images 
radiographiques, le modele calcule etant obtenu par 
deformation du modele generique de fagon que ledit modele 
calcule suive une forme la plus proche possible d'une 
isometrie du modele generique tout en maintenant en 
coincidence les reperes du modele generique deforme avec 
les reperes de controle reconstruits a l'etape (c) , 
caracterise en ce qu'au au cours de l'etape (a), les deux 
images radiographiques sont prises simultanement, par 
balayage, en deplagant en synchronisme, dans une meme 
direction de translation (3a) non parallele aux directions 
de prises de vues, au moins une source radioactive (5, 11) 
emettant deux faisceaux de rayons ionisants (10, 16) 
respectivement dans les deux directions de prise de vue (7, 
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13) . 

2. Procede selon la revendication 1, dans lequel : 
au cours de 1'etape (b) , certains des reperes de 

controle identifies, dits reperes de controle non stereo- 
5 correspondants ne sont visibles et identifies que sur une 
seuie image, 

et au cours de l T £tape (c), la position 
g6ometrique de chaque repere de controle non stereo- 
correspondant (C7-C25) dans le referentiel a trois 
10 dimensions est estimee a partir du modele generique, en 
deplagant les reperes de controle non stereo-correspondants 
du modele generique, chacun sur une droite joignant : 

d r une part, la source radioactive (5, 11) a 
l'origine de 1' image radiographique dans laquelle une 
15 projection de ce repere de controle non stereo- 
correspondant est visible et identifiable, 

et d' autre part, ladite projection de ce repere sur 
1' image radiographique, 

les reperes de controle non stereo-correspondants (C7-C25) 
20 etant ainsi deplaces jusqu'a des positions respectives qui 
minimisent la deformation globale du modele generique de 
l'objet a observer, 

3. Procede selon la revendication 2, dans lequel au 
cours de l'etape (c) , on minimise la valeur de la somme 

25 quadratique : 

S = AXM* ( -*,o) 2 , 

ou X est un coefficient constant, m est un nombre entier 
de ressorts fictifs reliant chaque repere (C1-C25) du 
modele generique a d'autres reperes de ce modele, ki est 
30 une valeur de raideur predeterminee du ressort fictif 
d'indice i, x io est la longueur du ressort fictif d'indice 
i dans le modele generique initial, et Xi est la longueur 
de ressort fictif d'indice i dans le modele generique 
de forme. 

35 4. Procede selon l'une quelconque des 
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revendications precedentes, dans lequel : 

- au cours de l'etape (b), au moins certains des 
reperes de contr61e identif iables sont des reperes de 
controle stereo-correspondants (C1-C6) visibles et 
identifies sur les deux images, 

- et au cours de l'etape (c) , la position 
geometrique des reperes de controle stereo-correspondants 
(C1-C6) est directement calculSe a partir de mesures de 
position des projections de ces reperes sur les deux 
images . 

5. Procede selon la revendication 1, dans lequel au 
cours de l'etape (b) , on repere sur chaque image 
radio graph i que des lignes de contour correspondant a des 
limites de l'objet observe et/ou a des lignes de plus 
grande densite optique a l'interieur desdites limites, ces 
lignes de contour comprenant des projections des reperes de 
controle sur les images radiographiques . 

6. Procede selon la revendication 5, dans lequel au 
cours de l'etape (c) , on determine des reperes du modele 
generique correspondant aux reperes de controle, lesdits 
reperes du modele generique comprenant des portions dudit 
modele generique qui se presentent tangentiellement par 
rapport aux rayons issus des sources radioactives et ayant 
generd les images radiographiques. 

7. Procede selon la revendication 6, dans lequel 
l'etape (c) comporte les sous-etapes suivantes : 

(cl) crder un modele generique recale en adaptant 
la taille du modele generique et la position de ce modele 
generique dans le referentiel, pour que les projections 
respectives du modele generique recale a partir des deux 
sources radioactives correspondent sensiblement aux deux 
images radiographiques, 

(c2) selectionner des reperes du modele generique 
dont les projections sur au moins une des images 
radiographiques a partir de la source radioactive 
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correspondante, sont les plus proches des lignes de contour 
reperees au cours de l'Stape (b) , 

(c3) d6finir une surface enveloppe formee par des 
rayons issus de chaque source radioactive (5, 11) et ayant 
contribue a generer lesdites lignes de contour des images 
radiographiques, 

(c4) determiner certains reperes du modele 
generique recale correspondant a des surfaces dudit modele 
generique recall, qui sont tangentes auxdites surfaces 
enveloppes, les reperes du modele generique recale ainsi 
determines correspondant aux reperes de controle, 

(c5) et determiner la position geometrique de 
chaque repere de controle par projection du repere 
correspondant du modele generique recale, sur la surface 
enveloppe correspondante. 

8. Procede . selon I'une quelconque des 
revendications precedentes, dans lequel les deux directions 
de prise de vue (7, 13) sont perpendiculaires l'une a 
1 1 autre. 

9. Procede selon la revendication 8, dans lequel on 
fait emettre par chacune des sources radioactives (5, 11) 
un faisceau de rayonnements ionisants (10, 16) dans un plan 
perpendiculaire a la direction de translation (3a) . 

10. Procede selon l f une quelconque des 
revendications precedentes, dans lequel les deux faisceaux 
de rayons ionisants (10, 16) sont emis respectivement par 
deux sources radioactives (5, 11) . 

11. Dispositif d'imagerie radiographie pour la raise 
en oeuvre d ? un procede selon l'une quelconque des 
revendications precedentes, ce dispositif comportant : 

des moyens d 1 emission de rayons ionisants 
comprenant au moins une source radioactive (5, 11), ces 
moyens d 1 emission etant adaptes pour emettre respectivement 
deux faisceaux de rayonnements ionisants (10, 16) vers un 
champ d'observation (4) contenant au moins un objet (20) a 
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observer, dans deux directions de prise de vue (7, 13) non 
parall£les, lesdits moyens d f emission etant deplagables 
selon une direction de translation (3a) non parallele aux 
directions de prise de vue de fagon que les deux faisceaux 
5 (10, 16) balaient simultan£ment le champ d f observation, 

au moins deux dispositifs de detection (6, 12) 
disposes respectivement face aux deux faisceaux de rayons 
ionisants (10, 16), au-dela du champ d T observation (4), 
pour mesurer les rayonnements ionisants ayant traverse 

10 ledit champ d' observation, ces deux dispositifs de 
detection etant d6plagables en synchronisme avec les moyens 
d' emission dans ladite direction de translation, 

des moyens (17) pour prendre simultanement au 
moins deux images radiographiques a deux dimensions du 

15 champ d' observation, par balayage simultane du champ 
d f observation avec les moyens d 1 emission (5, 11) et les 
detecteurs (6, 12) dans la direction de translation (3a), 

des moyens d' identification (17) pour 
identifier, sur chaque image radiographique, des reperes de 

20 contr61e predetermines (C1-C25) appartenant audit objet a 
observer, 

des premiers moyens de reconstruction (17) pour 
determiner une position geometrique de chaque repere de 
controle dans un referentiel a trois dimensions (X, Y, Z) , 

25 a partir d'un modele generique predetermine correspondant 
audit objet, ce modele generique comportant des reperes qui 
correspondent aux reperes de controle identifies sur les 
images radiographiques, 

et des seconds moyens de reconstruction (17) 

30 pour calculer la forme a trois dimensions d'un modele 
representant ledit objet a partir du modele generique, 
lesdits seconds moyens de reconstruction etant adaptes pour 
determiner le module calcule par deformation du modele 
generique de fagon que ledit modele calcule suive une forme 

35 la plus proche possible d'une isometrie du modele generique 
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tout en maintenant en coincidence les reperes (C1-C25) du 
modele generique deforme avec les reperes de controle 
reconstruits par les premiers moyens de reconstruction. 

12. Dispositif selon la revendication 11, dans 
5 lequel les moyens d* identification sont adaptes pour 
identifier des reperes de controle (C7-C25) visibles et 
identif iables sur une seule image radiographique, dits 
reperes de controle non stereo-correspondants, et les 
premiers moyens de reconstruction sont adaptes pour estimer 

10 la position geometrique des reperes de controle non stereo- 
correspondants en dfeplagant les reperes non stereo- 
correspondants (C7-C25) du modele generique, chacun sur une 
droite joignant : 

d'une part, la source radioactive (5, 6) a 

15 l'origine de 1' image radiographique ou une projection de ce 
repere de controle non stereo-correspondant est visible et 
identifiable, 

et d' autre part, ladite projection de ce repere sur 
1' image radiographique, 

20 les premiers moyens de reconstruction etant adaptes pour 
deplacer ainsi les reperes de controle non st§reo- 
correspondants jusqu'a des positions respectives qui 
minimisent la deformation globale du modele generique de 
l'objet a observer. 

25 13. Dispositif selon la revendication 11 ou la 

revendication 12, dans lequel les moyens de rep6rage sont 
adaptes pour reperer sur les deux images certains des 
reperes de controle (C1-C6) qui sont visibles et 
identif iables sur lesdites deux images, dits "reperes de 

30 controle stereo-correspondants" , et les premiers moyens de 
reconstruction sont adaptes pour determiner la position 
geometrique des reperes de controle stereo-correspondants 
par calcul a partir de mesures de position des projections 
de ces reperes sur les deux images. 

35 14. Dispositif selon . l'une quelconque des 
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revendications 11 k 13, dans lequel chaque detecteur (6, 
12) comprend une ligne (6a, 12a) de cellules de detection 
perpendiculaire a la direction de translation (3a), les 
faisceaux de rayonnements ionisants (10, 16) etant 
5 perpendiculaires a ladite direction de translation. 

15, Dispositif selon l'une quelconque des 
revendications 11 a 14, dans lequel les moyens d' emission 
(5, 11) et les detecteurs (6, 12) sont deplagables sur une 
distance de balayage d f au moins 70 cm. 
10 16. Dispositif selon l'une quelconque des 

revendications 11 a 15, dans lequel les moyens d' emission 
comprennent deux sources radioactives (5, 11) k l'origine 
respect ivement des deux faisceaux de rayons ionisants (10, 
11) . 
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3D reconstruction method using non stereo-corresponding 2D-contours on X-rays: 
case of the femur 

S. Laporte, D. Mitton, A. Gabrion, E. Brugere, J.B. Roy, W. Skalli, Laboratoire de Biomecanique, Service de 
Chirurgie orthopedique (Paris, Amiens Nord, F) 



Introduction: Existing 3D reconstruction methods for spine and pelvis using stereo-radiography are 
now quite accurate, but they require identification of the projection of specific anatomic points on the 
image (1 ). In the case of bone structure like the knee, there are not enough anatomical points to 
apply these methods. The aim of this study was to develop a 3D reconstruction algorithm based on 
identifiable non stereo-corresponding 2D-contours on the X-rays. The feasibility of the method was 
assessed in the specific case of the distal femur. 

Methods: Calibration of the X-rays environment and determination of the localization of the X-ray 
sources is a prerequisite, as well as a rough estimation of the object shape and location. Therefore 
the gross object is defined using 3D coordinates of N points in the X-rays environment referential. 
The method to get the reconstructed object from the gross one can be broken up in 5 steps: 

- For each view (Vi) and image (li), cones (Lij) are constructed based on each identified 2D-contour 
(Cij) and on the source X-rays location (Si). 

- From the approximate object location, the N 3D-points are projected on li (let Pi be a projected 
point). Then 3D points are considered as potentially generating a given 2D-contour Cij if the distance 
from Pi to Cij is lower than a given threshold. Subset of P'ij 3D points are therefore defined for each 
view. 

- From P'ij points subset, we extract P"ij points that correspond to a change of curvature on the 3D 
object, and therefore really yield a bone 2D-contour Cij on the image. 

- Those P M ij points (pertaining to the gross object) are then projected on their associated cones (Lij) , 
yielding 3D-points M"ij that pertain to the reconstructed object. 

- Finally , the M"ij 3D points are used as control points to get the N points of the reconstructed object 
by a kriging technique. 

For preliminary evaluation, 3 knees were stereo X-rayed in a calibrated environment (lateral and 
anteroposterior incidences) at 'CHU Amiens Nord*. One stereo-corresponding and 3 non-stereo- 
corresponding points were used to get the gross object shape and location, from a generic femur 
(556 points). The reconstructed objects were then compared with CT-Scan 3D reconstructions of 
these femurs (1 mm slices). 

Results: Feasibility of the method was assessed. Mean error for each distal femur was about 2.5 
mm, showing that a realistic 3D model can be obtained even with very few punctual anatomic 
information. 

Discussion: 3D reconstruction of bones is faced with the difficulty of 2D contours which do not 
correspond in the 2 views, and do not correspond to a predefined anatomic 3D contour. The 
proposed method allowed to get an acceptable solution. Work is in progress to improve the accuracy 
and to apply this method to other bones. 

1 Mitton et al MBEC Vol38(8) p1 33-1 39 
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A method for determining a three -dimensional form of a 
body from two-dimensional projection images includes pro- 
ducing at least one projection image of the body and 
producing corresponding model projection images from a 
generic model of the body. The body projection images and 
the model can be input into a computer- assisted processing 
unit and the generic model can be adjusted to the body by 
comparing image contents of the body projection images 
with the model projections images at a two-dimensional 
level only. The three-dimensional form of the body can be 
determined from the adjusted generic model. 
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METHOD AND DEVICE FOR DETERMINING A 
THREE-DIMENSIONAL FORM OF A BODY FROM 
TWO-DIMENSIONAL PROJECTION IMAGES 

RELATED APPLICATION DATA 

[0001] This application claims priority of U.S. Provisional 
Application No. 60/492,715, filed on Aug. 5, 2003, which is 
incorporated herein by reference in its entirely. 

FIELD OF THE INVENTION 

[0002] The present invention relates generally to the field 
of image production for anatomical structures in order to 
visually support medical and/or surgical operations. It finds 
particular application in conjunction with determining a 
three-dimensional form of a body or anatomical structure 
from two-dimensional projection images. 

BACKGROUND OF THE INVENTION 

[0003] Computer-assisted, stereotactic systems, which 
operate with the aid of body structure data obtained from 
tomographic detection systems and with the assistance of 
x-ray images produced in situ, are known, for example, from 
U.S. Pat. Nos. 4,791,934 and 5,799,055. Supporting an 
operation using x-ray image recordings is discussed, for 
example, in U.S. Pat. Nos. 5,967,982, 5,772,594 and 5,784, 
431. 

[0004] Even if accurate medical navigation is provided, 
the current prior art still operates with the aid of body 
structure data originating, for example, from tomographic 
detection systems, such as computer tomography devices or 
nuclear spin tomography devices. The patient to be treated 
is thus positionally registered in situ with respect to the 
previously ascertained image data. Operating instruments 
are then virtually displayed in the same relation to the image 
data as to the actual patient, in order to make the body 
structure data or possibly also x-ray image data useful to the 
surgeon in the operating theater. 

[0005] One drawback of such methods, in which tomo- 
graphic recordings (CT, MR) or x-ray images are produced 
especially for navigation within the context of a treatment, 
is the radiation load on the patient. Another drawback is the 
high cost, since such devices are very expensive both to 
purchase and to maintain and operate. 

[0006] Attempts have been made to develop systems that 
may be used without body structure data caplured separately 
in advance, for example, on the basis of statistical models of 
image data sets for body structure. However, such systems 
lack the required accuracy for the patient to be treated in 
each case. 

[0007] DE 100 37 491 Al and WO 99/59106 describe 
methods for providing 3-D information with the aid of 
fluoro-recordings. The starting point in all the methods is 
that of producing transillumination recordings of the patient, 
and/or the desired structure. A localization system is typi- 
cally used for this, in order to obtain spatial information with 
the recordings. DE 100 37 491 Al initially uses two fluoro- 
images, from which to reconstruct a rough 3-D model. 
Further recordings from different angles are used to specify 
the model even more exactly. In accordance with WO 
99/59106, at least three fluoro-recordings of the patient are 
generally made. ITiese can include anterior-posterior, lateral, 



and anterior- posterior with the head inclined back. In addi- 
tion to the transillumination recordings, photographs of the 
patient are also used. In this prior art method, the model is 
elaborately adjusted in three-dimensional space. 

[0008] A technique is known from WO 01/22368, using 
which a three-dimensional model is reconstructed using 
fluoroscopy images. The basis for this is a statistical model 
of the desired structure and at least one fluoroscopy image. 
The idea is to first position the model using the back- 
projected contours of the structure, such as are to be found 
in the fluoroscopy images. The model is then deformed and 
the contours of the deformed model are compared with the 
contours found in the fluoroscopy image. The quality of 
match is determined by calculating the error in the difference 
between the model contour and the fluoroscopy image 
contour. The two steps are repeatedly performed, i.e. the 
model is positioned and "morphed" until a good match is 
found. 

[0009] One drawback of this method is that the degree of 
match between the model and the fluoroscopy image data is 
determined in the three-dimensional space, which includes 
the model. Such calculations are very complicated and 
time-consuming, which compromises its practical applica- 
tion. Furthermore, accuracy suffers from the fact that this 
method only uses a few features of the fluoroscopy images, 
namely only the contour lines. 

SUMMARY OF THE INVENTION 

[0010] I n accordance with one aspect of the present inven- 
tion, the invention is directed to a method for determining a 
three-dimensional form of an anatomical body from two- 
dimensional projection images. The method can include 
producing at leas! one projection image of the body and 
producing corresponding model projection images from a 
generic model of the body. The body projection images and 
the model projection images from the generic body model 
can be input into a computer-assisted processing unit and the 
generic model can be adjusted to the body. The three- 
dimensional form of the body can be determined from the 
adjusted generic model. 

[0011] In accordance with another aspect of the invention, 
the invention is directed to a device for determining a 
three-dimensional form of a body from two-dimensional 
projection images. The device can include an image forming 
device, which produces at least one projection image of the 
body, and a computer-assisted data and/or image processing 
device, which produces model projection images corre- 
sponding to the body projection images from a generic- 
mode 1 of the body. The body projection images and the 
model projection images can be input into the data and/or 
image processing device using a data transfer means. The 
data and/or image processing device can include a compo- 
nent for adjusting the generic model to the body, wherein the 
form and, if necessary, the position of the generic model is 
adjusted, and wherein the match between the generic model 
and the body is determined at the two-dimensional level by 
comparing the image contents of the body projection images 
with the model projection images. The data and/or image 
processing device can also include a component for deter- 
mining the three-dimensional form of the body from the 
adjusted generic model. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0012] These and further features of the present invention 
will be apparent with reference to the following description 
and drawings, wherein: 

[0013] FIG. 1 is a schematic representation of a device for 
determining a three-dimensional form of an anatomical body 
from two-dimensional projection images in accordance with 
the invention; and 

[0014] FIG. 2 is an exemplary schematic diagram for a 
digitally reconstructed x-ray image used in conjunction with 
the present invention. 

DETAILED DESCRIPTION OF THE 
INVENTION 

[0015] In accordance with one aspect of the present inven- 
tion, the invention is directed to a method for determining a 
three-dimensional form of an anatomical body or structure 
from two-dimensional projection images. The method can 
include producing at least one projection image of the body 
and producing corresponding model projection images from 
a generic model of the body. The method can further include 
producing a generic model of the body and inputling the 
body projection images and the generic model into a com- 
puter-assisted processing unit. The generic model can be 
adjusted to the body, wherein the form and, if necessary, the 
position of the generic model can be adjusted. The match 
between the generic model and the body is only determined 
at the two-dimensional level, for example, by comparing the 
image contents of the body projection images with the 
model projection images. The three-dimensional form of the 
body is determined from the adjusted generic model. 

[0016] The match between the generic model and the body 
can be assessed at the two-dimensional level. When model 
projection images, which are two-dimensional in nature, are 
produced from a generic body model, these can be produced 
very quickly using appropriate hardware or software. In 
addition, comparing the model projection images to the 
likewise two-dimensional body projection images can be 
optimized for time and calculation efficiency. This can lead 
to a result far more quickly than with reconstructing and/or 
comparing matches in three-dimensional space. Further- 
more, less expensive hardware can be used. 

[0017] The body projection images can be x-ray images. 
The model projection images can be reconstructed projec- 
tion images, such as digitally reconstructed x-ray images 
(digitally reconstructed radiographs=DRRs), which origi- 
nate from volume or surface data sets for the model. The 
respective properties and advantages of using volume or 
surface data sets for the model are explained in even more 
detail in the following. 

[0018] In accordance with one embodiment of the inven- 
tion, comparing the image contents of the body projection 
images and the model projection images includes comparing 
outer contours of the projection images and/or features of 
the projection images within the outer contours. Comparing 
the image contents of the body projection images and the 
model projection images can include comparing depth infor- 
mation in the projection images. By using the image con- 
tents listed above when adjusting the generic model, it is 
possible to achieve a better match with the actual body 
and/or actual anatomical structure. 'Hie features in the inte- 



rior of the contour and/or the depth information contain 
additional criteria. The invention recognizes that projection 
images provide, not only data in one incision plane, such as, 
for example contour data, but also comprise additional 
information on the "depth," which can be used to optimize 
adjusting the generic model. Nonetheless, the image con- 
tents can be adjusted and compared at the two-dimensional 
level, which, in practice, allows a far quicker solution. 

[0019] In accordance with one embodiment of the inven- 
tion, the projection parameters of the body projection 
images are determined when producing the body projection 
images. This can be accomplished by means of a tracking 
system for the body and a recording device for the body 
projection images in order to produce the corresponding 
model projection images. Such tracking systems or naviga- 
tion systems can be (optically based) patient and device 
tracking systems, such as are understood by skilled artisans. 

[0020] In accordance with one embodiment, the generic 
model can be adjusted by way of statistical deformation. It 
can also be adjusted at least partially interactively with the 
aid of user inputs, wherein a free-form deformation model is 
used. 

[0021] For producing corresponding model projection 
images, a generic surface data set model can be used, whose 
data can then be used when adjusting the model. Alterna- 
tively, a generic volume data set model can be used. Also, 
surface data set models and volume data set models can be 
used as alternatives or together without departing from the 
scope of the present invention. 

[0022] According to another aspect of the invention, the 
invention is directed to a device for determining a three- 
dimensional form of an anatomical body from two-dimen- 
sional projection images. The device can include an image 
forming device and a computer-assisted data and/or image 
processing device, wherein the computer-assisted data and/ 
or image processing device can include components, such as 
hardware or software components, which perform one or 
more steps of the method described above. 

[0023] FIG. 1 shows a schematic representation of a 
device for determining a three-dimensional form of an 
anatomical body or structure from two-dimensional projec- 
tion images. 

[0024] In one embodiment, a (digital) x-ray recording can 
be taken of an actual body structure, for example, a vertebra 
1. A reference element 10 can be situated on the vertebra 1, 
using which the spatial position of the vertebra 1 can be 
detected by a tracking system and/or a medical navigation 
system. Such a navigation or tracking system is described 
more fully in co-owned U.S. Pat. No. 6,351,659, which is 
incorporated herein by reference in its entirety. The tracking 
system is schematically shown by a camera array 9 con- 
nected to a data or image processing device 7 via a data line 
8. The data or image processing device 7 is schematically 
shown here as a computer. A mobile or stationary image 
forming device, such as a C-arc x-ray device, is indicated by 
the reference numeral 2. ITie mobile or stationary image 
forming device 2 can include a radiation source 3 and, at the 
other end of the C-arc, an image intensifier 4, which likewise 
bears a reference element 5, which adds positional informa- 
tion to the recorded images. The image data can be trans- 
mitted to the computer 7 via the data line 6. 
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[0025] The computer 7 can include, among other things, 
image processing hardware/software, a CPU and data 
memories and a screen output, llie computer 7 can also store 
the data of the generic model, which is to be adjusted by 
means of the x-ray images produced by the x-ray device 2, 
in order to determine the ihree-dimensional form of the 
vertebra 1. 

[0026] For determining the three-dimensional form of a 
body structure in this way, digitally reconstructed x-ray 
images can be produced, which correspond to the x-ray 
images produced by the x-ray device 2. 'llie navigation 
and/or tracking system can ensure that the x-ray images of 
the vertebra 1 exhibit the same projection parameters (for 
example, location of the radiation source, location of the 
image intensifier, radiation parameters) as the digitally 
reconstructed x-ray images (DRRs). FIG. 2 schematically 
shows such a DRR being produced. Irradiating the vertebra 

12 from the generic model is simulated in the computer by 
means of a virtual radiation source 11. The model can be a 
surface model or a volume model. Using the radiation path, 
shown as dotted lines, a digitally reconstructed x-ray image 

13 can be created. The reference numeral 14 indicates a 
sample beam, which passes through the model. Using such 
a beam and/or using the sum of such beams, depth infor- 
mation in the DRR as well as information on features within 
the contours of the model can also be obtained. 

[0027] In accordance with one embodiment of the inven- 
tion, a three-dimensional form can be predicted and/or 
determined from the information of a plurality of two- 
dimensional projection images, for example, from at least 
one fluoroscopy image. In this way, a "model" of the 
observed form or structure (for example, a femur or a 
vertebra) is generated. This model can include statistical 
knowledge on variations in the form in a cross-section 
through a number of samples. However, this is not abso- 
lutely necessary. The form can be represented cither by a 
surface data set or, if more exact details of the inner structure 
arc also important, by a volume data set, or by a combination 
of both. In any event, a "generic form" of the structure of 
interest is provided. 

[0028] Changes to the form then can be coded in, for 
example, by (changing) a few form parameters. These form 
parameters can be, for example, the types of variation of a 
dot distribution model. They can also be node positions of 
free-form deformation boxes, which include the form. 

[0029] Furthermore, the images are at least roughly cali- 
brated. The projection parameters have to be known partially 
or entirely. This is enabled using image data, which may be 
determined in its spatial position via an image tracking 
system. 

[0030] Once the body projection images have been pro- 
duced and it has been possible to obtain corresponding 
model projection images from the generic model, the 
generic model is moved to a suitable starling position, 
whereupon the body projection images (e.g., fluoroscopic 
images) are compared with the model projection images 
(e.g., DRRs) rendered from the generic model. Ihe position 
of the model is, if necessary, corrected (i.e., the DRRs and 
the fluoroscopy images are assimilated to each other by 
rotating and shifting). A further step includes correcting the 
model form, wherein the DRRs and the fluoroscopy images 
are further assimilated to each other by changing one or 



more of the form parameters. The DRRs can be generated 
either software -based or hardware-based, if the generic 
model consists of a volume data set. If the model consists of 
a surface data set, then a DRR-like image can quickly be 
rendered, wherein the surface models are rendered as DRR 
images on standard graphic hardware. 

[0031] The position and form can be corrected for the 
model cither manually or automatically. Furthermore, they 
can be corrected either in sequence, in succession or in 
combination (simultaneously). 

[0032] One aspect of the invention also allows switching 
between surface model processing and volume model pro- 
cessing at any time during processing. For example, the first 
part of adjusting could be performed using the quicker 
surface model. Once a sufficient accuracy has been 
achieved, the result can be improved using the volume data. 
If an interactive mode is desired, the surface model would 
arguably be the better choice, and the resultant deformation 
can then be applied to the volume data. 

[0033] Deforming a statistical model purely intuitively is 
difficult to perform. To overcome this problem, an approach 
can be chosen in which a free-form deformation is used, in 
order to deform the model interactively or automatically. 
This deformation can be converted into statistical deforma- 
tion information, which allows switching between statisti- 
cally-based deformation and a free-form deformation at any 
time. One advantage here is that the user could prepare an 
initial deformation using intuitive free-form deformations, 
and the initial deformation can then be used as a starting 
point for statistical, automatic adjusting. Means can also be 
provided to improve the results of the automatic adjustment 
method, subsequently or at any time in the adjustment 
process. 

[0034] DRRs generated from the model are used and 
compared to fluoroscopy images. This results in better 
adjusting, since the depth information in the fluoroscopy 
images can also be used. If this additional information is 
ignored, more images may be needed, since the contours 
alone do not necessarily provide enough information 
(depending on the angle of the fluoroscopy images and the 
structure in question). 

[0035] The model can be assimilated to the fluoroscopy 
images automatically. This includes "morphing" the model 
and correctly positioning it. Landmarks can be detected or 
manually positioned, and the resulting model can be sub- 
jected to fine-tuning by using automatic registering and a 
morphing algorithm. It is also possible to only use DRR- 
fluoroscopy assimilation, in order to morph and position the 
model without the user intervening. 

[0036] There are DRR reconstruction methods that use a 
tetrahedron approach in order to model the volume data. The 
individual tetrahedrons are used to represent the various 
types of bone structures. The idea is to reduce the quantity 
of data necessary to display the volume while still providing 
the same information. This idea can also be implemented in 
accordance with the invention. 

[0037] Otherwise, the different bone densities can be 
approximated in the surface model, using different layers. 
One advantage of this is that everything can be morphed and 
rendered using conventional graphics cards, without special 
processing requirements. MTie additional layers, needed to 
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make the surface model more exact, can be automatically 
generated by segmenting additional iso-surfaces or also by 
scaling the model. 

[0038] The methodology described herein enables the 
form to be determined in accordance with practice, with a 
reasonable input of lime and calculating, and provides a 
better adjusted model than has hitherto been possible. 

[0039] It is to be appreciated that the method and the 
device described herein can include any of the features 
described therein, individually or in any combination. 

[0040] Although the invention has been shown and 
described with respect to a certain preferred embodiment or 
embodiments, it is obvious that equivalent alterations and 
modifications will occur to others skilled in the art upon the 
reading and understanding of this specification and the 
annexed drawings. In particular regard to the various func- 
tions performed by the above described elements (compo- 
nents, assemblies, devices, systems, compositions, etc.), the 
terms (including a reference to a "means") used to describe 
such elements are intended to correspond, unless otherwise 
indicated, to any element which performs the specified 
function of the described element (i.e., that is functionally 
equivalent), even though not structurally equivalent to the 
disclosed structure which performs the function in the herein 
illustrated exemplary embodiment or embodiments of the 
invention. In addition, while a particular feature of the 
invention may have been described above with respect to 
only one or more of several illustrated embodiments, such 
feature may be combined with one or more other features of 
the other embodiments, as may be desired and advantageous 
for any given or particular application. 

What is claimed is: 

1. A method for determining a three-dimensional form of 
an anatomical body from two-dimensional projection 
images, said method comprising: 

producing at least one projection image of the body; 

producing corresponding model projection images from a 
generic model of the body; 

inputting the body projection images and the model 
projection images from the generic body model into a 
computer-assisted processing unit; 

adjusting the generic model to the body; and 

determining the three-dimensional form of the body from 
the adjusted generic model. 

2. The method as set forth in claim 1, wherein adjusting 
the generic model to the body includes adjusting at least one 
of (i) the form and (ii) the position of the generic model. 

3. The method as set forth in claim 1, wherein adjusting 
the generic model to the body includes comparing image 
contents of the body projection images with the model 
projection images. 

4. The method as set forth in claim 3, wherein comparing 
image contents of the body projection images with the 
model projection images is determined at only a two- 
dimensional level. 

5. The method as set forth in claim 1, wherein the body 
projection images are x-ray images. 

6. The method as set forth in claim 5, wherein the model 
projection images are digitally reconstructed x-ray images 
which originate from volume or surface data sets for the 
model. 



7. The method as set forth in claim 3, wherein comparing 
the image contents of the body projection images and the 
model projection images includes comparing at least one of 
(i) outer contours of the projection images and (ii) features 
of the projection images within the outer contours. 

8. The method as set forth in claim 3, wherein comparing 
the image contents of the body projection images and the 
model projection images includes comparing depth infor- 
mation in the projection images. 

9. The method as set forth in claim 1 , wherein producing 
the model projection images includes determining projec- 
tion parameters of the body projection images while pro- 
ducing the body projection images. 

10. The method as set forth in claim 9, wherein the 
projection parameters are determined using a tracking sys- 
tem for the body and a recording device for the body 
projection images. 

11. The method as set forth in claim 1, wherein the generic 
model is automatically adjusted by way of statistical defor- 
mation. 

12. The method as set forth in claim 1, wherein the generic 
model is adjusted at least partially interactively with the aid 
of user inputs, wherein a free-form deformation model is 
used. 

13. The method as set forth in claim 1 , wherein a generic 
surface data set model is used for producing corresponding 
model projection images, data from the generic surface data 
set model being used when adjusting the model. 

14. The method as set forth in claim 1, wherein a generic 
volume data set model is used for producing corresponding 
model projection images, data from the generic volume data 
set model being used when adjusting the model. 

15. The method as set forth in claim 1, wherein a generic 
volume data set model and a generic surface data set model 
are used together or alternatively for producing correspond- 
ing model projection images, data from the generic volume 
data set model and the generic surface data set model being 
used when adjusting the model. 

16. A device for determining a three-dimensional form of 
a body from two-dimensional projection images, said device 
comprising: 

an image forming device which produces at least one 
projection image of the body; and 

a computer-assisted data and/or image processing device 
which produces model projection images correspond- 
ing to the body projection images from a generic model 
of the body, wherein the body projection images and 
the model projection images are inputted into the data 
and/or image processing device using a data transfer 
means, said data and/or image processing device 
including: 

a component for adjusting the generic model to the 
body, wherein the form and, if necessary, the position 
of the generic model is adjusted, and wherein the 
match between the generic model and the body is 
determined at the two-dimensional level by compar- 
ing the image contents of the body projection images 
with the model projection images; and 

a component for determining the three-dimensional 
form of the body from the adjusted generic model. 

17. The device as set forth in claim 16, wherein the image 
forming device is a digital x-ray device. 

***** 



